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EVALUATION 


The  effort  provided  an  initial  capability  to  create 
and  update  a  bathymetric  data  library  employing  voice  data 
entry  augmented  digitizing  and  editing  functions.  This 
effort  significantly  enhances  the  human  factors  aspect  of 
large  volume  geographic  point  data  digitizing.  The 
mathematical  data  manipulations  have  been  structured  to 
significantly  reduce  error  buildup  in  the  manipulated/ 
transformed  data.  Lastly,  the  digital  data  library  provides 
timely  query  responses  and  data  availability.  Future 
plans  are  to  increase  the  scope  of  the  mathematical  data 
manipulations  to  provide  for  a  wider  range  input  and  output 
formats,  to  introduce  a  wider  range  of  logical  retrievals 
at  the  data  base  for  increased  responsiveness  and  to  optimize 


le  voice  wtry  human  factors  aspects. 
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1  .  1  l'ui^'>2i** 

rllc  ot  tlllii  dlH  lunollt  IJJ  tWotoKl.  Kusl,  u  ji-oOUJlt.  1  L, 

l>ii>v  1  d«->.l  >.)f  ilio  subsyi>ti."m  and  lo<.tini>al  i •» i>li Ic-nis  which  were  encount¬ 

ered  in  l'ha;3e  1  ol  Bl>KS.  Second,  a  critical  account  and  evaluation  id  ‘Jiveii 
of  the  bolutions  which  were  provided  by  Synei  ticb  c‘ori>orat  ion  to  these  i>ioblenis. 

1 . 2  Subjects  to  be  Discussed 

This  document  is  divided  into  four  parts.  Part  I  is  dedicated  to  a 
discussion  of  problems  whose  solutions  are  mathematical  in  nature.  Part  11 
is  comprised  of  a  discussion  of  the  software  accomplishments  pertaining  to 
subsystem  problems.  Part  III  deals  with  the  myriad  contexts  in  which  the 
perennial  problems  of  accuracy  occur  in  data  input  and  manipulation.  Finally, 
Part  IV  critically  evaluates  the  scope  and  limitations  of  the  hardware  con¬ 
figuration  utilized  in  conjunction  with  the  system  software  provided  by  Data 
General  Corporation. 
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MA  TH  i;MA  r  U'Al,  Al>  I*  1 .  U  A  T  1 .  >N: ; 


j .  1  r^i 

tliiii  auction  la  to  ptovidi;  a  i.lotailuJ  account  ot  tliu 
maOiumat  leal  al^KU  1 1  lima  and  then'  ai>(>l  icat  ii>iu>  within  tlie  bUKS  inoblem  en¬ 
vironment.  The  thii'e  contexts  in  which  al*)oii  thins  have  been  developed  ale 
1  e>j  IS  1 1  a  t  ion  ,  I'ooidinate  t  vans  tonnat  ions  and  ijeoiiiaphn'  sectioniinj. 


Revjl  St  t  at  l<.>n 


2.J.1  Tlie  ptoblem  to  be  solved  by  n  point  re9istration  is  twofold, 
first,  a  function  must  be  constructed  whose  domain  is  comprised  of  physical 
X-Y  locations  on  a  dujitizin^  table  and  whose  ran9e  is  comprised  of  X-Y  cooid- 
inates  in  on  earth  iectan9ular  frame,  yiven  that  n  values  (3  to  8)  of  the 
function  are  known.  The  construction  of  this  function  is  referred  to  as  day 
*1*  n-point  re9i:;iration.  Second,  a  function  must  be  constructed  which  maps 
table  X-Y  locations  to  talile  X-Y  locations  9iven  that  n  (from  3  to  8)  values 
of  the  function  ,ire  known.  This  construction  is  entitled  day  'n*  n  point 
reqistration . 

Day  '1'  retjistration,  from  an  empirical  point  of  view,  results  in  the 
9eodetic  significance  of  points  on  a  map  which  has  been  placed  on  the  digit¬ 
izing  table.  Day  'n'  registration  results  in  mapping  back  to  day  1  the 
points  on  the  map  as  it  lies  on  the  table  translated  or  rotated  relative  to 
the  day  1  registration  of  the  map.  Since  both  algorithms  are  identical  from 
a  mathematical  point  of  view,  it  suffices  to  give  a  presentation  of  day  1 
registration. 

2.2.2  Method  of  Solution 

The  idea  of  the  n-point  registration  is  to  convert  X-Y  table  values  to 
a  form  which  can  be  converted  to  lat/long  values  on  the  earth's  surface.  A 
pure  table  X-Y  value  is  physically  meaningless. 

TTie  problem  is  a  typical  statistical  one  of  finding  a  function  given  that 
you  know  some  of  the  ordered  pairs  that  the  function  must  satisfy.  In  regis¬ 
tration  you  know: 

2-1 

^  y  10 


k 


(a)  tho  tablo  X-Y  valuis  of  the  registration  points 

(b)  the  map  scale  mil  values  of  the  lat/long  values  of  the  registration 

points  (the  map  scale  mil  values  are  found  by  taking  the  output 
from  the  map  projection  itself  and  converting  to  units  of  mils) 

To  be  more  precise,  let  X^,  be  the  map  scale  mil  X  values  of  the 

registration  points,  let  Y^,  Y',...y''  oe  the  map  scale  mil  Y  values  of  the 

registration  points,  let  X,,  X^,...X  be  the  known  table  X  values,  let 

1  2  n 

Y,,...Y  be  the  table  Y  values.  You  know  for  each  ordered  pair  of  X-Y  table 
In 

points  5(^,Y^  the  map  scale  mil  equivalent  X  ,Y  .  This  constitutes  the 
known  relations  which  must  be  approximated  by  the  function  being  sought.  What 
properties  should  the  function  have  beside  fitting  the  known  data?  Simplicity 
requires  that  the  function  be  linear.  Geometrical  considerations  suggest 
finding  two  function  f,q  such  that 


(a) 

f(X 

1  1 

(b) 

X  i 

Given  the  assumption  of  linearity  we  get 

(a)  f (X  ,Y  )  =  AX  +BY , +C  =  X^ 

11  11 

(b)  g(X  ,Y  )  =  DX  +EY.+F  '  Y^ 

11  i  1 

Now  we  must  impose  some  conditions  which  allow  us  both  to  solve  for  A,B,C,D,E 
and  F  unambiguously  and  which  lend  geodetic  meaning  to  the  physical  table 
X-Y  values.  That  is,  given  an  arbitrary  table  X-Y  pair,  the  functions  f  and 
g  must  map  the  pair  into  an  approximately  correct  map  scale  mil  ordered  pair. 

The  classical  approach  to  a  statistical  problem  in  a  linear  model  is  the 
least  squares  best  fit  algorithm.  The  idea  is  to  minimize  the  sum  of  the 
squares  of  the  differences  between  known  data  values  and  the  values  predicted 
staistically .  In  registration  we  have  known  map  scale  mil  values.  We  want 
the  functions  f  and  g  to  yield  predicted  map  scale  mil  values  such  that  we 
minimize  the  sums  of  the  squares  of  the  differences  between  the  predicted 
map  scale  mil  values  (i.e.,  f(X^,Y^),  g(X.,Y,)  and  the  known  map  scale  mil 
values  (i.e.,  X^,y^). 

More  explicitly,  given  that 

(a)  X^  =  AX.+BY,+C 

.11 

(b)  Y^  =  DX,+EY,+F 

1  1 
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(d  ,  (X^  -  (AX  +By  +01' 

1=  i  11 


II 

(d)  i;  iv^  =  (DX.  +EY.  +F)  J' 


i-1 


by  pickiiuj  the  correct  values  of  A,B,C,D,E,  and  f 


Let  H(X  =  5^  [X^-(AX^+By^+C) 


i=l 


i  ** 

and  G(y  .X  .y.)  =  z  Iy^-(AX.+By  +C)]^ 
^  ^  i=i  11' 


the  requirement  of  minimizing  these  e]q>ressions  is  logically  equivalent  to 
the  following  six  equations: 


(1) 

M  - 

3A  “ 

(2) 

M  - 

3B  “ 

(3) 

OH  _ 

3C  ~ 

(4) 

^  - 

3d  ~ 

(5) 

dG 

3e  “ 

(6) 

3G 

3f  “ 

Taking  the  derivatives  and  siiiplify4ng  we  get  the  following  six 
equations : 


n  2  n  n  n  . 

(1)  AS  (X  )^  +  bE  X  y  -  cE  X  -  E  x^x, 

i-l  1-1  ‘  ‘  i-1  i-i  ^ 


n  n  _  n  ”  < 

(2)  A  E  x.y  ♦  B  E  (y  )“'  +  c  E  y  -  E  x‘y 

i-1  ‘  ^  i-1  i-1  ^  i-1  ^ 


n  n  n  . 

(3)  A  E  X  +  B  E  y  +  n  -  E  x‘ 
i-1  i-1  ‘  i-1 


n  2  I'  n  n  . 

(4)  D  E  (X  )  +  E  E  X  y  ♦  r  E  X  -  E  y*x 

t-1  ^  i-1  ^  ‘  i-1  ‘  i-1  ^ 
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To  solvu  for  y 


the  invorso  of  X;  X 


n  1 

1 

1  ’  1 

1-' 

to  IJut 


(11  A  x'*;'. 

B 
C 


(2) 


U 


V 


Tho  oMttix  X  ^  oxiiitti  if  itx  detonxinAnt  ia  non-t«ro  which  sauMi  to  work 
alwdyx.  It  It  doex  nut  wuik,  new  reqlatretloit  points  must  bs  chosen. 

The  suftwAte  tut  this  Approach  requires  solvlnq  for  the  elements  of  X.S.W, 
tindinq  X  And  multiplyinq  X  ^12  And  X~^V.  HAVlnq  found  these  vsluas  ws  hsve 
tor  All  arbitiAry  tAhle  XV  pAlr  thAt 

(1)  X^  -  AX^tBY^+C 

li)  if  OX^*KY^tr 

Upon  conversion  of  x\y^  to  meters  «ra  multiply  by  the  mep  ecele  to  qet 
U^.vS  the  earth  scale  meter  equivalents  of  X^,Y^.  This  pair  U^.V^  Is  the 
input  into  the  Inverse  map  projection,  the  output  Is  the  pair  6,i|i  called 
the  latitude  and  lonqitude  of  the  table  point  X^«Y^. 


t  ;iv  111  1  > ,  ,  iiiii  1 1 1 :  1 1 1 1 Y  (1 1  /Viij  ii  i>i>i  1 1 

'lau  avki|'liil  .i|'|  iii.iili  lo  t>'i|  1  st  I  >11  ion  w>i:<  il«  vo  li  i  i  il  'iii-lion,i  '  jt  Syini  l  li  s 
Conniiiit  loll  111  1 1^1 1  t.  I V  II I  .11  ,  till'  ot  .■.olvili'l  fill  I  liu  X  ijllil  Y  il  I  ini  ii:,  1  out. 

I  llvic'|H,‘IKli‘lll  1  V  1  Ol.U  1  I  III  111  ini  II  1  III!  ..  i  l|i|  otlol;,  wlilill  l  DU  I  li  IC'liUlf  fiom  toifujlfj- 

ztiiij  siiihIo  iDiitlol  (loiiit:!  wlioii  a  reqiatiatioii  i iiii.i.  Cuiif  alili.* . 

J  .  J  .  4  Kiniii  rioal  Koaultii  and  Kvaluat  1 1 ii i 

SiiK'ti  tliu  matliomat  ica  I  modol  of  rugistiation  ia  boili  atatistical  and 
linoar,  tliero  aro  oiiiiu  i  i  oa  1  siLuationa  in  whioh  probloma  can  arise.  In  the 
ideal  situation,  tlie  source  map  is  not  distorted  badly  by  shrinkage  of  expan¬ 
sion.  We  liave  noted  tiiat  such  distortion  is  rarely  compensated  for  within  the 
linear  model.  A  distiuted  iiuip  can  be  registered  within  the  residual  tolerance 
level.  Tills  sliovild  not  be  taken  as  evidence  that  the  statistical  fit  has 
adequately  compens.ited  for  the  distortion.  The  model  is  linear)  theiefore, 
the  .scaling  which  occurs  is  uniform  along  an  axis.  Thus,  the  actual  map  shrin¬ 
kage  is  spread  luiiformly  over  the  map.  This  could  make  the  functions  con¬ 
structed  in  legist lation  both  meet  the  residual  test  and  fail  to  provide  an 
adequate  mapping.  Thus,  the  chart  must  be  known  to  be  in  good  condition  or 
poor  empirical  dat.i  from  the  table  might  bo  allowed  to  infiltrate  the  data 
base . 


It  is  crucial  Lo  pn'k  reg i st i at  ion  points  intelligently.  If  B  registra¬ 
tion  points  are  [>ii-ked  representative  of  the  map  as  a  whole)  that  is,  in  no 
obvious  geometrical  p.ittern,  and  the  document  is  in  good  physical  condition, 
then  the  results  of  registration  can  be  trusted.  All  mathematical  cximputat ions 
are  performed  in  over-kill  double  precision. 


2.2.5  Future  Topics  of  Interest 


Since  the  mathematical  modol  of  registration  is  linear,  the  idea  obviously 
arises  of  utilizing  a  model  in  which  non-linear  teims  cxintribute  to  statistical 
prediction.  The  linear  model  is  quite  satisfactory  given  a  chart  in  goixi  con¬ 
dition  and  an  intelligent  choice  of  reglsliation  points.  Thus  the  employment 
of  a  non-linear  model  would  bo  motivated  if  an  attempt  were  made  to  registei 
charts  in  less  than  optimal  condition. 


.CVS 


;s.i- 
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The  linear  model  could  also  be  supplemented  with  a  normally  distributed 
error  term  whose  expectation  value  is  zero  and  whose  standard  deviation  re¬ 
flects  user  input  in  accuracy  in  entering  a  control  point.  However,  it  seems 
arbitrary  to  assign  such  a  stamdard  deviation  for  an  arbitrary  user.  Thus, 
this  term  does  not  appear  in  the  functions  constructed  by  registration. 

2 . 1  Coordinate  Trans  format  ions 

^.3.1  Hrolileni  to  be  Solved 

The  problem  to  be  solved  by  the  coordinate  transf ormations  is  to  convert 
between  earth  rectangular  and  geodetic  coordinates  utilizing  the  Mercator  and 
Transverse  Mercator  projections.  The  conversion  from  geodetic  to  earth  rectan 
gular  coordinates  is  effected  by  the  map  projection  itself.  Inverses  had  to 
be  constructed  as  well  to  convert  earth  rectangular  to  geodetic  coordinates. 

2.3.2  Method  of  Solution 

The  projection  types  utilized  are  conformal  projections.  As  Thomas  has 
shown  in  "Conformal  Projections  in  Geodesy  and  Cartography",  Coast  and  Geodeti 
Survey  Special  Publication  251,  all  conformal  mappings  of  the  spheroid  upon  a 
plane  are  expressed  by  the  analytic  function 

X+iy  -  f(A_^ir)  where 
r  -  In  tan  | 

L  ‘l+5sln^  '  J 

and  ^  equals  the  latitude  of  the  point,  X  the  longitude  and  (  the  eccentrrciiy 
The  form  of  the  Mercator  and  Transverse  Mercator  mapping  is  then  determined  by 
which  line  or  lines  in  the  projection  are  to  be  held  true  to  scale  and  by  the 
necessary  geometric  form  of  map  elements  corresponding  to  meridians  and  paral¬ 
lels.  Once  the  form  of  f  is  determined,  the  real  and  imaginary  parts  of  Xtiy- 
f(X^ir)  are  equated  which  must  in  tprn  satisfy  the  wall  known  Cauchy  Riemann 
equations : 


(1) 

?x  - 

iX 

3X 

3r 

(2) 

3x  •> 

-  ^ 

ir 

3X 
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to  finally  obtain  x  as  a  function  of  X  and  r  and  y  as  a  function  of  X  and 
r.  The  mapping  equation  is  thus  derived. 

2. 3.2.1  Mercator  Mapping  Equation 


The  Mercator  projection  is  the  simplest  and  most  basic  of  the  conformal 
map  projections.  The  form  of  f  is  linear  and  the  initial  conditions  quite 
siRf>le.  As  Thomas  indicates «  the  initial  condition  is  that  the  scale  of  the 
projection  must  be  true  at  the  equator.  Thus,  for  a  latitude  of  i^O,  n«ln(l) 
■  1^,  and  Y“0;  and  x=aX,  the  product  of  the  equatorial  radius  and  the  longi¬ 
tude  X.  Thus  since  x+iy«f(X+ir)  we  have  aX+i. (^f (X+i. «f (X)-a (X+ir) .-x+iy . 
Equating  real  and  imaginary  parts  of  x+iy«a. (X+ir)  we  have 


(1)  x-aA 

/  f 

(2)  y“ar»  a.  In  ltan(7f  +  1(^1)  (l-^9in<f>) 

'42  l+4sin^'  I 


=  a.  ln(  (sin<>tl) 


(1-Csin<j>) 

<l+^siw$T 


The  scale  or  magnification  of  a  point  at  latitude  4*  given  a  conformal  pro¬ 
jection  is  the  Jacobian  of  x  and  y  with  respect  to  r  and  X  divided  by  Ncos^* 


where  N=< 


1  .  2, 

1-f,  sin  d 


the  distance  from  the  point  along  that  line  from  the 


minor  axis  of  the  spheroid  perpendicular  to  tl»o  line  tangent  to  the  point. 

In  this  case  the  scale  become  a  .  secdi. 

N 


2. 3.2.2  Inverse  Mercator  Mapping  Equation 

Given  the  Mercator  mapping  function  Y  =  fix)  which  maps  a  latitude  to 
an  earth  scale  meter  value,  find  a  method  to  calculate  the  latitude  given 
the  earth  scale  meter  value. 

The  Mercator  mapping  equation  in  question  ta)tes  tlie  form: 


<''/sCALAc’ 


[ZEl^C  /  2  I 
tan  t  m)  .  l-7.KCC.siN(r) 

4  2  1+ZECC.SIN(P)  J 


Where  Y  =  the  earth  scale  meter  value  and  P  =  the  latitude  value. 

Given  this  form  for  f (x) ,  an  expression  i,(X)  must  be  derived  such  that  the 
form  p  *  6(P)  is  obtained.  Given  an  initial  approximation  P„  to  P,  Weustetn 
iteration  is  used  to  imi'rove  the  desired  latitude  value. 


2-H 


To  derive  the  form  P  -  6(P)  we  proceed  as  follows.  Since  SCAIAC  is  a 
constant  we  eliminate  it  to  get 

ZECC/2  ) 


-  A. In 


IPl) 

2 


I-/.I  <  I  .SIN(l  ) 
1  t/lA  1-.  blN(l') 


where 


Uui 


Wan  (TT  + 

[  ^ 

-  A.  In  ^an  +  IPljj  +  Y(P) 

/  /.11  1'/ 2  \ 
Y(P)  -  A. In  (  '  1 

'  I  i/.iiCi  .;;1N(1-)  / 

./i  .  c.  1/2  III  /  1  -/.!■:>  t  .MN(P)  \ 

\  1 ♦ZtCf.SIN(P)  J 


Y  ( 1' )  -  A . 


Su  if  wtr  let  b  -  '/.t> 'i' .  b  I N  (P)  and  use  the  series  expansion  1/2  In  1*B 

1-B 


b*l/JB  ♦  J/‘jb  ♦...  with  AE2  =>  A.2ECC,  AE4  -  A.ZECC.ZECC/ 
AEb  *  A.Ztcr.ZtlC.ZEc'C/^j  we  yet 
Yd')-  S1N(P).(AE2  *  SIN^(P).(AE4  *  SIN^(P).AE6) 


31 


Therefore : 


Yj  “  A. In 

^tan 

p 

Yj  -  Y(P) 

In 

A 

p 

EXP  - 

Y(P) 

C  " 

♦ 

P  «  2  . /arctan 

y(p) 


Thus  we  have  derived  the  correct  form  P  -  6(P)  for  use  by  the  iteration  scheme. 

The  initial  approximation  P,  to  P  is  obtained  by  considering  ZECC  «  0i 
that  is,  a  spherical  earth  model  is  employed.  Proceeding  as  before  we  havei 


«  A. In 


P.  -  2. 


|arctan  ^  exp(yj/^)j  ~ 


Thus  P,  is  calculated  exactly. 

Having  obtained  P.  and  P  -  6(P)  Wegstain  iteration  is  used  as 

.  ;i  OUALlXi  i’KACIHiAiilfi 

follows:  -  ....  i  _ _  „ 


IWJ  oorx  t'UiXlilSiUii)  10  DOC 


\  \  ' 


(1)  Initializatiun 


P  =  P 
*^0  0 

=•  G(PJ 

?!  =  G(P^)  +  G(P^)  -  P^ 

"■  ■  r- 

'•'■V 

(2)  Al  1  li-i  4it  lull  ii  1 1 
Pm  1  -  (.;(Pn) 

Pm  1  -  i.i(Pml)  t  G(Pn-»-l)  -  Pn*l 
Fn»l  -  pn 
G(Pn+l)  -  Pntl 

(3)  Cuiivci ijuiicc  at  tolerance  E 
Piul  -  Pn 

- V  L 

Pnf  1 

Tliiti  will  occur  unleiiS  G^<P)  equals  1  at  P.  If  convergence  does  not 
occur,  ttie  initial  apinuximation  P„  is  output  as  the  latitude. 
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2. 3. 3. 3  Transverse  Mercator  Mapping  Equation 

As  Thomas  points  out,  the  initial  condition  of  the  Trwsverse  Mercator 
projection  is  that  the  scale  is  true  along  the  central  meridian  of  the  map. 
Thus,  for  a  point  at  longitude  zero  and  latliuile  4>  wc  have  x=U.  Given  that 

xtiy-f(A*ir)  w«?  Ii.ivo  iy-f(it)-i.  i the;  di  tit, nice  along  the  mer  idian 


t'rum  the  oguaior  at  longitude  \  to  the  point  at  latitude  4>  with  K=  the  radius 

j  2  2  -3/2  • 

of  curvature  ).a.(l-f,  sin  ♦)  '  .  but  r  -  \  R  sec4>d(|>  by  the  condition 

^  N 

on  r^or  all  confoimal  projections.  Tlius  dr-K  secitidi|>  and  Rd4=Ncos4>dr .  Tlierefore 
Ncosi^dr-f  (r)  sincel  is  zero.  Thus  we  have  x+iy=S4>=f  (r)  .  Now  we  want  to 
subject  X  and  y  to  the  Cauchy  Rieaiann  equations  so  that  equations  for  x  and  y 
in  terms  of  A  and  r  are  needed.  Thomas  proceeds  to  expand  xtiy^f (X'tir)  about 
the  point  z~ir  in  a  Taylor  Series.  Then  real  and  imaginary  components  are 
equated  yielding  the  desired  fonns  for  x  and  y  in  terms  of  X  and  r.  Once  all 
terms  are  computed  and  the  Cauchy  Riemann  equations  applied,  Thomas  reaches 
his  desired  mapping  equations. 


Tlio  scale  K  for  the  projection  at  a  point  X,^  is 


Ncos^ 


=  1  Xx  .  (1  +  tan^r)^'^^ 

Ncos<^  9X 

2. 3. 2.4  Inverse  Transverse  Mercator  Equations 

The  following  algorithm,  taken  from  Rapp  and  Sprinsky,  Page  30  is 


eB^>loyed . 


sec4  X 


(1  +  2t^^  ♦  ♦  1  X 

12  0  N 

a  2  2  2 

•♦■24j^  ♦  bhj  +  8tj^  Oj  ) 


4-4* _ _!  (1  +0.^)  X  ^  1  i  (1  *  ^^1^ 

2  ^  N  2  4 


(5  ♦  28t, 


2  ..  4 

♦  -4nj^ 


+  X  ^  _  (61  ♦  90tj^  ♦  *^*^1*  *  1070^^ 


7  2  0 


-  162e'‘ 


sin^4'  ♦ 


46j'^t  *sin%’)  X® 
^  n 
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where:  tj^  »  tan 

»  e'^  cos 

e'  »  second  eccentricity  squared 
-  tout  i<oint  Idtitudti 
X  UTH  Kabtinq 

I 

N  •  a'/ti  -e*  sAn^iJ)^ 
a*  -  .<J496.a 

a  >•  aemi -major  axis  of  the  ellipsoid 

A>  -  ioiiqitudu  di f t'oreiicu  from  cential  mutidian  of  point 
-  latitude-  ut  (joint 
e^  ->  eccentricity  squared 

Havinq  computed  AX,  and  the  longitude  itself  must  be  computed  as 
fol lows : 

X  -  -(AX  -  (6(J0  -  n)  -  3)  p  ) 

where:  X  ^  longitude 

p  -  .01745327252 
n  >  zone  number 

The  first  negative  converts  longitude  positive  eastware  to  positive 
westward. 

A  first  approximation  to  foot  point  latitude  is  computed  from: 

♦•-S(l  +8*  (AA  ♦  S*  (BB  ♦  (CC  ♦  S^DD)))) 
where : 

AA  -  -3e^/6 

BB  -  (12e^  ♦  45e'')/120 

CC  -  -(48e*  ♦  1023#"  ♦1170e‘)/5040 

DD  -  (192e'  *  18384e''  ^750996*  ^6048a' 1/362880 
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Y  -  Y/(.9996a(l-c^) ) 

Y  -  UTM  Nort)iiny 


.«| 'I 'I  tik  I  hu  l  I  i.>l  4*'  <*>>-'  ubt  j  I  lu.-d  by  Nc-wtuii  1 1  cl  a  I  lull , 


li>  : 


Ji*  -  -  - 

^(1*1)  ^  tl)  V" 


wild  Ml  K  -  Si  -  i>  J 

diiJ:  s>  ^  •  Aifi'  •  (b/^)  Hill  *  (C/4)  Min  (4<>') 

-  (U/b)  Mill  (bi^‘) 

A  -  ‘  '  b4‘’'  *  2S'(>^ 


“  ■  4‘’*  *  ib^  *  5U^ 

•  a--  • 

“  -  ih"' 


aiiil  I 


(1  -  M*  sin^^')^ 
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J  .  i .  i  : m  1  .  t  at iii  li  I  .  I » ' i tit  1 1  i»tt  (t 

'I'llf  'rilvim-ia  .It  I  iVt.l  I '  I  1 1  )  t  t  t  1 1  >li:t  llaVt-  lifiii  i  1  tjiit  txi;.  )  y  lialtil  ulitl  .ij  luVi-tl 

111  many  t  ai  t  tit|  i  ,t|  .it  i  ,  |<iti|iti;..  'I'lif  iiivii-.t-  it.  llif  Mtit  jitu  Mtt|.|  iitt|  wn-o  iIl'- 
vclojifti  .tt  byintlit:,  t  .  .1 1  .til  at  1  mi .  Tlit  luvt  i  :.t  itj  llif  I'l  .insvii  i.t-  Miitaloi  inaj.- 
piiitj  tia;.  liofii  iilitaiiittl  timii  Ka('|>  iiitl  Ltpi  1  iii.ky . 

J.  t.4  bunt  1 1  l  t.a  \  kfbii  1 1  aixt^  l:v.iluat  imi 

Tlio  1 1  alia  t  til  iiu  t  ioim  ftoni  ealtli  lectaiiqular  to  y«.-otji  al'lll  c  coonliiiates  , 
wlioii  comliinod  witti  tlu>  traiiafotniatioiit;  doii veil  in  registration,  ate  aecurate 
to  within  a  second  of  arc  wittiin  tlie  well-known  bounds  of  application  of  tlie 
map  projections  ttiemselves.  They  are  accurate  enough  to  satisfy  any  require¬ 
ment  which  stays  wittiiti  this  limitation.  The  inverse  Transverse  Mercator  map¬ 
ping  IS  confined  to  one  zone  at  a  time.  The  inverse  Mercator  mapping  equation 
handles  aibitrary  earth  rectangular  values. 

.  J .  5  Future  Topics  of  Interest 

The  inverse  transverse  Mercator  mapping  is  restricted  to  transforming 
earth  rectangular  coordinates  wittiin  a  single  zone  to  geographies.  It  might 
prove  beneficial  to  generalize  the  mapping  to  allow  zone  overlap  if  Uiis  is 
possible . 

2 . 4  Geographic  Sectioning  Algorithms 

2.4.1  Problem  to  be  Solved 

The  sectioning  algorithms  provide  the  capabilities  to  perform  circle, 
path,  and  polygon  searches  within  areas  generated  on  the  surface  of  the  earth. 
The  earth  is  imagined  to  be  a  sphere  of  unit  radius. 

2.4.2  Method  of  Solution 

To  perform  a  given  search,  the  routine  must  be  called  twicej  the  first 
call  generates  the  appropriate  parameters  to  describe  the  region  in  question, 
and  the  second  call  accesses  points  and  calculates  whether  they  fall  inside 
the  test  region . 

■  rViJi:  lo  hSJT  .IVAlliY  rRACIlCAHL^ 

■  >...i  Jil  i  i- uiOilijiii.J  10  ooc  . 
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.  I  K-  :it  ui  ‘  ji 

(A)  t  tit  I'.il  I 

i>li  til  l!  lull  lliii  lit|iut  lut-l(-)iii|  vuliicu  ut  tlie  cc-iitL'i  dlid  I'uiJlial 
I'OlIlt  die  eutivel  led  tu  a|ilici  I  cal  iwol  dl  Itate -a  oit  a  unit  ta^diere . 

(B)  Sucund  Call 

Kirdt  the  input  tedt  point  IN  convut tud  to  Nphuricdl  coordinatNN 
uit  d  unit  Nphete.  I'tieii  the  qruat  elide  dietaneeN  fru«  the  center  point  to 
the  radial  point  and  from  the  center  point  to  the  teat  point  ere  celculeted. 
If  the  latter  distance  ie  leee  then  the  former,  the  test  point  is  within  Che 
re9ion  defined  by  the  center  end  rediel  points. 

The  distances  ere  celculeted  es  followei 


Let  K  end  S  be  vectors  eminetiny  from  the  ori9in  of  the  ephere  to 
the  2  points  in  question.  The  dot  product  of  R  end  S  is  the  cosine  of  the 
en9le  between  Che  two  vectors.  The  en9le  is  the  9rset  circle  dietence  in 
question,  rormelly  we  heve 


R‘S 


but 


Xj  ■  sinS^'Cosf^ 
Y^  -  sind^'sin#^ 


•  cos6^ 


Xj  -  einOj  cosfj 
Yj  -  eina^  coefj 
Zj  «  cos6^ 


where  9  is  latitude  and  9  is  longitude.  Thus  R'S  *  sind^  sinR^  (coeia^-a^)) 
♦  cos9^  cosd^, letting  W  •  dietence  sought  we  hew* 


Tins  FA  SI" 


li.  I  ■;> 


\;;i  n  i 

XO 
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W  -  ARCCOS 


sinU^  siiiO,^ 


(cos  )  +  cosOj^  cosO^ 


Q.E.D. 


J  .  4  .  J  .  J  >  1  y  ^|on  Sea  t  ch 

(A)  b'llst  Ciil  1 

Un  titat  coil,  tlic  input  Idt-lonij  values  of  tlm  poIy9onal  vertices 
are  converted  to  Caitesian  X-Y-Z  coordinates  on  the  surface  of  a  unit  sphere 
centered  at  the  uri<jin  of  the  coordinate  system. 

(b)  Second  Call 


Un  second  call  the  input  lat-lon9  values  of  the  test  points  are 
converted  to  Cartesian  X-Y-Z  coordinates  on  the  surface  of  a  unit  sphere 
centered  at  the  ori9in  of  the  coordinate  system. 

To  calculate  whether  a  9iven  test  point  lies  within  the  polygonal 
test  region,  the  following  insight  is  utilized.  Suppose  the  n  vertices  of 


the  polygon  have  coordinates 


<X  ,Y  ,Z  *  and 
n  n  n 


the  test  point  has  coordinates<X^,Y^,Z^>.  Consider  the  planes  defined  by 
the  triangles  whose  vertices  are  ^  , 

<x,,»,.z,>} . 

<X^, Y^,Z^^andthe  line  defined  as  that  line  connecting  the  origin  to  the 


coordinates<X^, Y^,Z^>. 


Let  W^.Wj, 


W  he  the  points  of  intersection 
n-2 


between  the  line  and  each  of  the  planes.  Let  triangle 

and  a  given  plane,  the  angles  described  from  the  point  of  intersection 
to  the  triangular  vertices  ' 

If  3  ~  then  the  point  of  intersection  lies  within  the  triangular 

region.  If  this  occurs  for  at  least  one  triangle,  the  point  is  obviously 
in  the  polygonal  region. 
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The  following  is  a  formal  derivation  for  a  test  point  and  a  given 
triangle.  First  the  equation  of  the  plane  of  the  triangle  is  calculated,  then 
the  point  of  intersection  is  calculated,  then  the  origin  of  the  coordinate 
system  is  translated  to  this  point  of  intersection,  and  finally  the  angles  are 
calculated  and  summed. 

(1)  e'a  1 1  u  1  .i  t  lull  lit  plane  ol  liiungle 

'I'lie  ge'iiei.il  equation  lit  a  plane-  is  given  Liy  the  formula  A-X  +  b-Vt 
• '-li  *  L>- •  i'o  v.le- 1 1- 1  mi  lie-  tin:  planar  e-guatiun  given  1  points  ,Zj^  > , 

"Xj.Y^.Z,  ,  We-  have 


+  (X^Z^  -  Z^Xj) 


D  =  - 


^(Y^Z^-Z^Y^)  -  ^i<>‘2^3"^2^3^  ^  3''‘ 2*2 


Given  that  the  three  points  are  triangular  vertices,  the  above  coefficients 
determine  the  plane  of  the  triangle. 

(2)  Point  of  Intersection 

Suppose  that  RP(I),  1-1,2, 3  is  the  XYZ  coordinate  of  the  test  point. 
Then  the  equation  of  the  line  in  space  through  the  origin  and  the  test  point  is 

X  =  _Z_ 

RP(l)  RP12)  RP(3) 

Suppose  RP(1)  y  U.  Then  we  proceed  as  follows: 


y  =  RP(2) • X 
RP(l) 


Z  =  RP (3) .X 
RP(1) 

C-  [fy  (3)  -X 
RP(1) 


'  A-XtB-rKP(2)  -X  I  +  C-  jiy  (3)  -x"!  +  D 

[^RP(l)  J  I^RPd)  J 

►  X-  A  t  B.kP(2)  ♦  C-  RP(3) I  =  -  D 

L  RP(1)  RP(1)J 


THIS  PAGE  IS  BEST  QUALITY  PRAP.TTpAH^ 
COPY  FURBISHED  TO  DDC  _ 
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X 


=  _ 

'A-t^B'RP(2)  +  C-  RP(3) 

RP(l) 


i.ii  i;  n-Ki'i.') 

ki'  ( I ) 

K  -  l';  Kl'  (  I) 

Kl  ' (  I  ) 

H  -  am;ii 

tllOll  Wl'  IliiVl.' 

X  -  ;11 
M 

l.i-'t  KPl(l),  1-1,2,  i  b«:  the  coordinate  ot  the  point  of  intersection. 

Thun  we  have 

RPI  (1)  -  X  -  ^ 

H 

RPI(2)  -  RP(2)  -HPKl) 

RP(l) 

RPl(J)  =  KiM3)  -RPld) 

RP(l) 

(3)  Translate  origin  to  point  of  intersection 

Suppose  '■X^,Y^,Z^>  is  an  arbitrary  point  in  X-Y-Z  space.  Then 
SJK(I),  I.>1,2,3  is  defined  as  follows 

SJI(l)  =  X  -RPI(l) 

1 

SJI (2)  =  Y  -RPI (2) 

3 

SJI (3)  -  Z^-RPI (3) 

Then  SJI(I)  is  the  translated  point 
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A 


\ 


CB  =>  , 

A  -  RPI 

Note  the  CB  is 

the  array  A  translated  to  a  new  i 

column  of  CB  is 

a  coordinate.  To  find0  ,6  ,  and 

the  following  vector  dot  products 

(1) 

<CB(I.l)  .  CB(I.2)> 

(2) 

<CB(I,1),  CB(I,3)> 

(1) 

vCB(I,2) ,  CB(I,3)> 

according  to  the  formula 

^ » 

-t/  2  2  2 

lL2  „2 

A'B  ^ 

^  \/a  ♦  a  ♦  a 

UK  *  K 

V  X  y  * 

V  X  y 

whe  re 

A-A  C+A  + 

x-1  yj 

A  ? 

z*k 

BxB  C"  * 

B  C 

x.i  yj 

z*k 

Given  that  the 

following  assignments  are  made 

(1) 

NHAGl  ••  (CB(1,1)^  * 

CB(2,1)^  ♦ 

(2) 

RMAG2  -  (CB(l,2)^  + 

CB<2,2)*  ♦ 

(3) 

RMAG3  -  (CB(1,2)^  ♦ 

CB(2,3)*  ♦ 

z 


cos  ( 


,2,1/2 

,2,1/2 

.2.X/2 


we  have 


(1) 

CB(I.l) 

■  CB(I,2) 

-  RHACl  •  RMAG2  • 

cosBj  ■ 

DOTl 

(2) 

P’ 

CB  ( 1 , 1 ) 

’  CB(1,3) 

-  RMAGl  •  RMAC3  • 

cosB^  ■ 

DOT2 

13) 

P 

CB  ( 1 , 2 ) 

♦ 

•  CB(I,3J 

-  RHAG2  >  BMAG3  • 

cosO^  - 

DOT3 

Thus  each 
follows.  Take 
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Therofore 


(1)  10^  “  DOTl/(RMAGl -1^102) 

<2)  lOj^  “  DOT2/(RMAGl-RMAG3) 

(3)  103  =  [X)T3/(RMAG2 -RMAGa) 

Therefore 

(1)  -  AHt.GOii(lb^) 

(2)  *  AKCCUii(ia^) 

(3)  6^  ‘  ARCCOS(ia^» 

2.4. 2.3  Rdth  Search 
(A)  First  Cdll 

First  the  3  points  input  which  dsfins  ths  psth  srs  converted  to 
XYZ  coordinates  on  the  surface  of  a  unit  sphere  yielding 
<X^,Yj,Z^>.  The  vertices  of  the  quadrangle  defined  by  the  path  ere  found 
as  follows.  Let  AX-X^-Xj,  AY-Yj-Y^  and  Z-Z^-Z^.  The  four  vertices  ate 


(1) 

<X3.Yj.Z3> 

(2) 

<X2-AX,Y2-AY, 

%-AZ> 

(31 

<X3-AX,Y3-Ay, 

(4) 

<X3+AX,  Y^+AY 

,  Z3*AZ> 

Thie  follows  by  the  obvious  sysssetry  of  the  sphere 
<B)  Second  Call 

This  is  just  a  4  sided  polygon  search. 
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2 . 4 .  i  Suat  t-i-  uiui  llistury  of  Afproacii 


Tile  soLtioiiiny  aisjuritliiiis  were*  coinijletoly  develoiieii  in-liuubu  dt  Synectics 
Corporation. 


2.4.4  Enn>iricdl  Results  and  Evaluations 

The  sectioniny  algorithma  have  bet:n  tested  thoroughly  and  are  accurate 
in  typical  application  environments.  The  only  source  of  approximation  consists 
in  assuming  that  the  earth  is  circular.  Thus  geodetic  paths  on  the  ellipsoid 
are  considered  to  be  great  circle  paths  on  the  surface  of  the  earth.  For 
large  earth  test  areas,  this  could  result  in  some  inaccuracy  in  polygon  searches. 
The  extent  of  the  inaccuracy  has  not  been  rigorously  analyzed. 


2.4.5  Future  Topics  of  Interest 

The  polygon  sectioning  routines  are  employed  within  th^  context  of  de¬ 
termining  whether  a  given  quadrangular  region  overlaps  a  given  polygonal  re¬ 
gion.  An  algorithm  needs  to  be  developed  which  determines  whether  two  such 
regions  overlap.  At  present,  a  serial  search  method  is  employed  to  solve  this 
problem  and  this  seems  to  be  a  needless  consumption  of  time. 

It  might  also  be  of  interest  to  develop  an  algorithm  which  would  appro¬ 
ximate  the  region  of  overlap  in  question  with  a  many  sided  polygon.  In  this 
way,  a  third  region  could  be  introduced  so  that  a  determination  could  be  made 
whether  three  regions  overlap.  Thus,  in  general,  an  iterative  procedure  would 
exist  to  determine  when  n  regions  overlapped  and  an  approximate  polygonal  re¬ 
presentation  of  the  overlap  would  be  obtained. 

Finally,  the  polygon  search  might  be  generalized  to  more  complex  polygonal 
regions;  that  is,  regions  where  angles  could  be  obtuse  between  adjacent  sides. 
This  would  probably  be  required  if  the  iterative  procedure  described  above 
were  developed. 
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.';i  i-iiiiN  111 


;>i  'l■•l'WAKl:  Al  l  I  ini‘i.  1  iiIimi  n  is 


1.  1 

Itu.’  (Ull  I'l'L.c  ol  tills  .'■I'v.l  lull  Is  ti>  ills.'Uiis  stiiiu-  ut  tlic  mull’  ]  II 1 1’ I  cs  t  1  Ilsi 
ilovolupmuiil  s  wliii'li  luivu  oi  1  111  icii  liiii  iii.|  till’  iii.iliui,  ul  ilui  I'liasu  1  BIJKS  sott- 
ware.  Tlio  iJis,iiss  lulls  i  iiu  1  iiilu.i  aiu  nut  iiiti-inkil  lu  i  xliaiist  all  taui’ts  ut  the 
^y^tem.  Ilu’  tui'us  ut  I'vii'li  sepal  ate  liisi.'Ussiuii  will  be  tiuin  the  key  novel  eapd" 
bilitii’s  pioviiteJ  liy  the  suttwaie  in  satistyini)  the  le’iui  i  emeiit  s  of  the  thlee 
BDKS  subsystems;  di  ij  i  1 1 2a  t  i  on  and  voiee  entry,  bateh  and  data  base. 


i  •  ~  Uiiji  t  izdt  ion  and  Voiie  Entry  Subsystem 
3 . .  I  i'ilL'-* 

The  key  aoi  ompl  i shment  of  this  BDKS  Subsystem  is  that  of  provtdinij  the 
eapability  to  diijitize  soundinij  data  from  nautical  charts  and  store  this  sound- 
iny  information  in  a  compact  and  conveniently  accessible  form.  Soundiny  data 
may  be  entered  manually  with  a  keyboard  and  vocally  via  communication  between 
the  Eclipse  C300  and  the  Threshold  bOO  Voice  Kecoynition  terminal/display 
device.  As  the  voice  data  is  enteied,  a  visual  display  is  provided  on  the 
Tektronix  4010  display  unit. 

The  soundiny  data  may  also  be  edited  by  the  user  in  a  very  simple  way 
usiny  either  the  special  keyboard  or  the  'voice  box’.  Thus,  a  complete  capa¬ 
bility  is  provided  to  create  and  edit  features  of  sounding  data  which  makes 
maximum  use  of  system  resources  to  ease  the  task  of  the  user  in  data  entry  and 
edi  t . 


i.i.2 


Fat hoy rams 

An  important  capability  of  this  BDKS  subsystem  is  that  of  diyitiziny 
fathograms;  that  is,  graphs  of  depth/time  coordinates  supplemented  with  such 
parameters  as  geographic  fixes,  loxodrome  bearings,  and  ship  velocities.  Once 
a  digitized  fathoyram  file  has  been  created,  a  file  of  geographic-depth  posi¬ 
tions  can  in  principle  be  constructed  and  entered  into  a  data  base.  Thus, 
the  analog  infoim.it  ion  of  the  fathoyram  is  inteyratable  into  the  general  BDKS 
data  base  framwork. 
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Tins  rAi'.K  IS  hbST  QUALITY  rRAOHOMi^i 
.vuM  COi  Y  FUlvNlSUliD  TO  DDC  _ - - 


J.J.  )  K>-»|  I  ,-.1  i  .1 1  tun 

Till'  uliii’i  1 1  lull  I'l  ii-'l  I  Ml  I  at  lull  {luM  bc-i-ii  t  Ik<i i>ui|lil  y  lii  mi'IImmli-I  lli  t.c'i,ljijii 
J.J  ut  Hum  lii  m  luiii  II I  .  11  mIiuuIiI  liuwi-Vi-r  be  elil|>li.iM  i  ^eil  licli;  tliat  Lilt.'  UeVeluji- 

nKMit  aiii.1  linpK  iiii-iit  at  lull  ut  tlilM  a  lijoi  i  tlini  witliiii  tins  subsystem  pruvides  tile 
subsystem  user  a  veiy  i- 1 ;  i  i.  i  eii  t  and  easily  eini-loyed  metliml  ut  t  eiji  s  teri  tnj 
eharts  so  that  inti>rmat  ion  in  a  fuatuie  diijitized  fuim  that  ehail  is  easily 
accessed  and  aecuialely  edited.  In  the  Lineal  Input  System  (LIS)  no  such 
method  was  emp loyed land  the  errors  inherent  in  both  the  LIS  day  '1'  and  day 
'll'  reijist  ration  pioiedures  have- been  eliminated.  Within  the  conceptual  struct¬ 
ure  of  the  reijistrat ion  scheme  developed  by  Synectics  Corporation,  it  becomes 
an  easy  matter  to  edit  table  files  even  when  they  are  created  by  the  batch 
program  which  converts  arbitrary  geographic  files  to  table  files.  Thus,  an 
arbitrary  geographic  file  constructed  from  information  in  the  data  base  can 
in  principle  be  converted  to  table  form,  edited,  and  reconverted  to  geographic 
coordinates.  A  more  thorough  discussion  of  this  will  be  found  in  the  discus¬ 
sion  of  the  batch  program  which  converts  geographic  files  to  table  files,  in 
Section  3.3.1  of  this  document. 

3.2.4  Review  Mode 

The  capability  exists  to  display  the  data  of  a  table  file  on  a  Tektronix 
4010  display  given  a  user  selected  window.  Each  feature,  whether  consisting 
of  trace  data,  discrete  points,  or  soundings,  is  searched  and  data  which  falls 
within  the  window  is  displayed.  This  provides  the  user  with  all  he  needs  to 
discover  where  he  desires  to  make  edit  changes.  The  scope  of  such  iui  editing 
capability  in  relation  to  the  Data  Base  is  further  discussed  in  Section  3.4 
of  this  document. 

3. 3  Batch  Processes  Subsystem 

3.3.1  Geographic  to  Table  Conversion 

The  novel  feature  provided  by  this  capability  is  that  the  resultant  table 
file  can  be  immediately  edited,  in  principle,  by  employing  a  day  'n'  registra¬ 
tion.  This  can  be  seen  as  follows.  When  the  file  is  converted,  the  registra¬ 
tion  points  of  the  file  are  also  converted.  The  resultant  file  has  been  day 
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\  registered  in  the  sense  U\at  a  i^erfeot  fit  Is  obtained  with  least  square 
beat  fit  coefficients  of  >1>.  This  means  that  any  table  point  in  tlie  file  need 
only  be  scaled  to  earth  scale  meters  followed  by  a  translation  from  the  chart 
earth  scale  nk'ter  lv>wer  loft  hand  corner  to  a  true  oiiqin  to  obtain  the  true 
uaith  scale  meli't  value  ot  the  table  point.  The  plotted  I'hail,  when  placed  on 
the  table,  is  retilsteied  as  a  day  ii  reqistiallun  of  the  peifoi't  tit  file. 
Tlierefote,  qiveii  that  a  set  of  registration  i'OiiiLs  ate  in  the  qeoqrapliic  file, 
an  arbitrary  >jev>'iraphit.'  file  can  be  brought  to  table,  edited,  and  reconverted 
to  geographic  t<.iim. 

The  mathemat icai  simulation  of  this  is  as  follows.  Suppose  are  the 

latitude  and  longitude  of  an  arbitrary  coordinate  in  the  geographic  file,  bet 
X,  ,y,  be  the  map  scale  meter  value  of  the  lower  left  hand  corner  of  the  chart. 

la  Is 

Let  X.,y.  be  the  map  scale  meter  value  of  the  coordinate  (4,i).  The  table 

A  >P 

X,Y  of  the  point  is 

(1)  X  ■((X  -X  )  *  J  inchesi  .  conversion  to  mils 

A  la 

(2)  y  “Ky^-yj^)  ♦  3  inchesi  .  conversion  to  mils 

Let  X  ,Y  be  the  map  point  X.Y  (4«i)  when  the  map  is  placed  on  tlie  table 
n  n 

which  is  produced  by  a  plot  of  the  table  file.  Then  regietretion  will  map 

X  y  to  x,y  and  X,,y^  need  only  be  calculated  from  eguatlone  1  and  2  to  give 
n  n  '  ♦ 

the  map  scale  meter  values  of  When  scaled  to  earth  scale  meters  by 

multiplying  by  the  map  scale,  the  earth  scale  meter  value  le  obtained  for 

(#.X). 

3.  J..’  Table  to  Geographic  Convereion 

A  novel  feature  provided  by  this  capability  is  that  two  distinct  kinds  ot 
table  file  can,  in  principle,  be  converted  to  geographic  formj  digitised  table 
files,  aj)d  table  files  created  by  geographic  to  table  convereion.  The  dif¬ 
ference  between  the  conversion  le  that  the  letter  table  files  ere  transformed 
to  earth  scale  router  coordinates  as  described  in  Section  3.1.1,  whereas  the 
former  table  files  are  converted  to  eerth  scale  meter  coordinates  by  accessing 
the  day  I  least  square  best  fit  coefficients  created  when  tha  table  file  is 
legisterud  as  dusciibed  in  Section  3.3.3. 
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A  .•(vvtiul  IK'Vi'  I  ^it  ttiku  \‘vt|k  Vcl  k*  i  K>li  Ilk  \lu  W.t'^  kit  whkvil 

l>Oklkt  to  iMklll  vtut  .k  k  :•  tkv'.ltk'tl  kll  t  I  .ii  o  lu.tl  kkk  >:• .  roiiil  to  |t'liil  tial.t  t  >>ii- 
kllkltsk  V>t  tWt>  tkt'k.ti  1  kit  o  oiiklkotlkioii  (HkliitM  IK'V  III  t  lli>|  skii't'i'Ski  I  vit  ly  III  llio  tcatkiio. 

A  utiaikilii  1  kiu-  IS  I  Hip  I  1  o  1 1  1  y  asiiukiitsl  t>>  oxist  l>i>lwt'i.'ii  t  ho  two  points.  It 
the  data  is  >.\MiVk'k  ( ist  tk>  ktk.'s>.|i  aplitos  aii«i  thou  li.ivk  tk.>  tahlo  using  a  di  t  tL’iciiit 
m>ip  pio  ik'k't  iv>n,  I  lio  sliatght  lino  will  ivivoi  klittoiont  map  sk.ali>  motet  points. 
Tike  solkitit-Hk  to  till:!  pi'olilum  i.'oikiiiMts  Ilk  goiioiating  X-Y  points  on  thu  linu 
k'onneo  t  i  ng  ttio  twk.>  iKiints  in  the  kiriginal  tablo  tilu  at  inotements  ot  IJ?  mils 
in  oithoi  X  oi  Y.  As  an  i not omviit al  X-Y  |H>sition  tin  tliti  line  is  obtained,  it 
is  conveitod  tv>  govigi apliios .  Tims,  the  implioit  line  as  a  whole  is  lepicsont- 
atively  oonvettod  tki  vieogtaphios  in  the  talile  ti>  geinjiaphii-  oonvetsikgi  tunotion. 

i.  J.  J  Plot  Kkinot  ions 

Thu  plot  sottwaie  has  several  features  wvirthy  of  especial  notti'u.  The 
develo|Mnent  of  this  software  is  such  tha^  in  principle,  any  other  type  ot 
plotter  can  bti  added  to  the  BURS  hardware  conf igurat ion  sini>ly  by  accessing 
an  intermevliate  plot  file.  Second,  the  Xynetics  interfav.'e  routines  are  written 
in  simi>le  Fortran  and  are  in  principle,  usable  by  other  i.\M«|>uter  systems.  Fin¬ 
ally,  the  development  of  further  plot  functions  is  readily  integratable  into 
the  current  sviftware  due  to  the  single  calling  seviuence  of  any  plot  function. 

Fur  exan^le,  functions  existing  on  other  coim'uter  systems  interfacing  with  a 
CalcoRip  plotter  may  be  integrated  into  the  BOKS  system  to  produce  the  same  type 
of  output. 

J • <  Pata  Base  Subsystem 

J.4.1  Sect  ickiiing 

The  Data  Base  subsystem  provides  the  capability  of  creating  BDRS  geogra¬ 
phic  files  on  the  basis  of  search  criteria  which  are  thoroughly  discussed  in 
Section  i.4.2  of  this  document.  In  principle,  the  following  sequence  of  steps 
could  be  taken,  utilising  this  capability,  to  quickly  edit  the  files  included 
in  the  source  sectioned  files.  First,  the  created  geographic  file  is  converted 
to  table  and  plotted.  Hie  plot  is  then  reuistered  and  the  table  file  reviewed 
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until  llu  il.il.i  l>’  1«‘  I'viilfil  IS  tuUnil.  'I'lu'  Lilil  iiiiidi'  Is  tliLii  k'titcit.'iJ  tu  tdil 
the  t  ill'.  I  .n  It  li'.iluie  islitisl  nieiitltie:  the  biiuiee  111  anti  tioe uiiien t  nuiiiber 
of  the  Stiuiie  SIS  I  loiievl  tile,  as  well  a:i  the  number  nt  the  feature  relative 
to  this  stiuitu  Sit.  l  1  eiietl  tile.  The  edited  f,le  is  Linn  lei'on  Vi:  I  ted  to  qeo- 
qraphies  with  a  tla>|  set  at  each  tiatuie  whi.h  has  been  edited.  t'inally,  the 
source  settioin’d  files  are  updati'd  to  refleet  the  edits  that  were  made  on 
the  table  tile.  The  folKiwiiuj  diaqram  reflects  the  flow  of  this  sequence. 

SOUKCE  SECTlONtD  TILKS 


GEtXlRAPHIC  FILE 


TABLE  FILE 


□  □  □  □  □ 

1 

□ 

1 

□ 

I 

□  □  □  □  □ 

EDITED  SOURCE  SECTIONED  FILES 


EDITED  TABLE  FILE 


GEOGRAPHIC  FILE 


-c-T  , 
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TllK  I'KOBLtMS  OK  ACCUKACV 


4. 1  Pur|)use 

The  KuriJode  of  this  section  is  to  anelyze  the  critical  contexts  within 
which  considerations  of  accuracy  are  relevant.  In  each  such  context,  the  in¬ 
tent  is  to  identify  both  where  considerations  of  accuracy  have  been  analyzed 
and  where  further  analysis  is  required. 


4, 2  Data 

4.2.1  Accuracy  of  Source  Analog  Data 

4.2. 1. 1  Charts 

The  accuracy  of  data  portrayed  on  man-made  navigational  charts  and  maps 
is  clearly  suspect  in  many  instances.  Charts  become  distorted  because  of 
shrinkage  and  expansion,  and  the  distortion  is  clearly  a  highly  complex  non¬ 
linear  phenomena  from  a  mathematical  point  of  view.  The  algorithm  of  registra¬ 
tion  does  little, if  anything,  constructive  to  con4)ensate  for  this  distortion 
as  registration  is  represented  within  the  classical  linear  model  of  maximum 
likelihood  which  assumes  uniform  stretching  or  shrinking  in  an  arbitrary 
direction. 

When  such  charts  are  digitized  and  converted  to  geographic  coordinates, 
information  becomes  available  to  the  data  base  which  seriously  threatens  its 
integrity.  Thus,  the  ability  to  edit  the  data  base  becomes  of  paramount  im¬ 
portance.  A  serious  look  should  be  taken  at  the  model,  provided  by  Synectics 
Corporation,  in  Section  3.4.1  of  this  document,  which  purports  to  encompass 
such  an  editing  capability. 

The  accuracy  of  data  portrayed  on  charts  produced  mechanically  by  plotters 
on  such  systems  as  the  Lineal  Input  System,  presents  emother  problem  as  the 
giaterial  out  of  which  the  charts  are  made  is  much  more  resistaint  to  distortion. 
These  problems  will  be  discussed  momentarily  in  Section  4.3  of  this  document. 


li 
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<1 . 2 .  1 . 2  Data  Entry  of  Analog  Source  Materials 

Data  is  cMitored  by  tlu;  user  within  the  digitization  subsystem  of  BDKS. 
When  graphical  feature  data  is  entered  with  a  cursor,  the  data  is  accepted  as 
is  by  the  system;  that  is,  there  is  no  error  term  emi^loyed  by  the  software  to 
smooth  out  the  errors  which  occur  in  data  entry  itself.  In  a  linear  model, 
the  typical  assumption  of  statisticians  is  that  s'  'h  error  can  be  modeled  by 
a  probability  distribution  with  a  zero  expectation  and  a  variance  characteris¬ 
tic  of  the  particular  user.  Whether  the  inclusion  of  such  an  error  term  would 
be  fruitful  within  BDRS  itself  has  not  been  analyzed  empirically. 

A  natural  question  to  ask  is  how  to  check  whether  the  data  entered  by  the 
user  is  representative  of  the  data  supplied  on  the  chart.  The  proof  plot  pro¬ 
vides  the  check  required.  If  the  proof  plot  is  a  perfect  overlay  of  the 
original  source  chart  as  it  is  taped  on  the  table  when  digitized,  then  each 
digitized  point  is  within  an  epsilon  of  its  correct  position.  That  is,  de¬ 
pending  upon  the  resolution  of  the  data  when  digitized  (from  1  mil  to  10  mils) 
and  accuracy  of  the  plotting  device.  It  may  appear  (to  the  human  eye)  that 
two  chart  points  overlap  each  other,  but  in  reality  may  actually  be  a  few 
mils  apart. 

The  accuracy  of  the  sounding  data  entered  by  the  user  can  be  checked  by 
the  user  as  he  enters  the  sounding  value.  This  value  is  displayed  on  the 
Tektronix  4010  display.  If  the  user  is  employing  the  Threshold  500  Voice  unit 
his  entry  is  further  displayed  on  the  Threshold  display  unit. 

4.2.2  Accuracy  of  Processed  Data 

4. 2. 2.1  Registration 

The  algorithm  of  registration  is  thoroughly  discussed  in  Section  2.2.2 
of  this  document.  In  this  section  several  corollaries  pertaining  to  the  pro¬ 
blem  of  accuracy  will  be  deduced  from  the  mathematical  properties  of  the 
algorithm. 

It  was  noted  in  Section  2.2.4  of  this  document  that  registration  points 
should  be  picked  judiciously;  that  it,  representative  of  the  chirt  as  a  whole 
and  not  in  any  obvious  geometrical  pattern.  The  reason  for  this  latter  con¬ 
dition  finds  its  justification  in  the  mathematical  fact  that  an  earth  scale 
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Thus,  it  IS  the  users  i eS))ons ibi 1 i ty  to  ensure  that  he  is  not  allowing  poor 
data  into  the  BUHri  system.  There  is  a  check  which  can  be  employed  for  this 
case.  It  is  discussed  in  Seclion  4 .  .1 . 4: .  2  of  this  document. 

A  final  coiollary  of  this  algorithm  is  that  the  proper  registration  of 
a  good  source  chart  geneiates  an  exiieinely  accurate  mapping  from  the  coordinate 
frame  of  the  table  to  the  earth  rectangular  frame  in  day  1  registration.  The 
simplicity  of  the  algocitlim  conjoined  with  the  par V.ti oiling  of  syiratietric 
matrices  to  find  invetses  results  in  a  very  accurate  mapping. 


4 . 2 . 2 . 2  Coordinate  Transformations 

When  geographic  coordinates  are  mapped  to  the  earth  scale  meter  frame  and 
then  back  to  the  geographic  frame,  the  resultant  geographic  coordinates  are 
within  a  second  of  arc  of  the  original  geographic  coordinates.  The  case  of 
iterating  this  procedure  to  check  for  ai,  accumulation  of  harmful  error  has 
not  been  checked  rigorously. 


Since  the  coordinate  transformations  and  the  registration  transformation 
are  so  accurate,  the  following  sequence  constitutes  an  excellent  test  of  whether 
the  source  analog  chart  is  distorted  badly. 


(1) 

(2) 

(3) 

(4) 


Register  the  chart  and  build  a  table  file. 
Convert  table  file  to  geographic  file. 
Convert  geographic  file  to  table  file. 

Plot  table  file. 


^iv  A 
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If  the  plot  overlays  clecuily  with  the  source  data,  then  the  chart  Is  in 
good  condition.  Otherwise,  the  accuracy  of  the  chart  is  seriously  suspect. 

4. 3  Problem  of  Validation 

In  order  to  evaluate  the  accuracy  of  the  map  projections,  charts  must 
be  utilized  which  represent  data  to  the  degree  of  accuracy  being  tested. 

The  Hydrographic  Center  has  provided  such  a  chart  for  the  Mercator  projection 
test.  Charts  of  such  high  caliber  will  also  be  provided  for  each  mapping 
demanded. 

To  evaluate  the  accuracy  of  the  polygon  search  algorithm,  a  high  pre¬ 
cision  Gnomic  projection  chart  is  required.  The  great  circle  vertices  of 
the  polygon  map  to  straight  lines  on  such  a  chart.  Thus,  one  may  construct 
with  precision  arbitrary  earth  polygonal  regions  on  the  chart  and  test  arbi¬ 
trary  points  for  inclusion  within  the  defined  chart  polygonal  region. 

If  the  accuracy  of  the  map  projections  is  precisely  known,  the  accuracy 
of  the  functions  constructed  in  registration  can  be  tested  by  simply  com¬ 
paring  their  output  with  the  known  output  of  geographic  points  either 
identical  to  or  different  from  the  registration  points. 

The  overall  problem  of  all  of  these  validation  procedures  is  that  they 
demand  charts  which  are  'known'  to  be  accurate  within  the  degree  of  pre¬ 
cision  for  which  a  test  is  to  be  made.  If  the  charts  are  machine  produced 
from  the  Lineal  Input  System  (LIS) ,  then  the  projections  must  conform  to 
those  employed  in  LIS  itself.  Whether  this  is  desirable  or  not  has  not 
been  rigorously  analyzed  within  the  BDRS  research. 

The  procedure  to  this  point  has  been  to  implement  the  LIS  geographic 
to  table  projections  for  the  Mercator  and  Transverse  Mercator  mappings. 

The  equations  have  been  checked  and  validated  within  BDRS.  Whether  such 
a  procedure  will  be  effective  for  dealing  with  the  Polyconic  and  Leunbert 
Conformal  projections  is  not  at  present  known. 
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SECTION  V 


SYSTEM  CONFIGURATION 


5. 1  Purpose 

The  purpose  of  this  section  Is  three-fold.  First,  to  provide  a 
detailed  description  of  the  hardware  confiijuration  employed  in  Phase 
I  of  the  BDRS.  Secondly,  to  identify  ^tnd  discuss  key  technical  areas 
which  pertain  to  the  development  of  the  BDRS  functional  capabilities. 
Thirdly,  to  present  conclusions  and  recommendations  pertaining  to  the 
BDRS. 

5.2  Hardware  Configuration 

Figure  No.  5-1  illustrates  the  BOC  hardware  configuration  which 
consists  of  the  following: 

o  Data  General  ECLIPSE  C300  Processor  -  128K  core  memory 
o  Data  General  Magnetic  Tape  Units,  9  Track  (2) 
o  Data  General  6012  CRT 
o  Centronics  Line  Printer 

o  Data  General,  92MB  Disk  Drive,  dual  portable 
o  Data  Automation  Digitizing  Tables  (2) 
o  Tektronix  Graphic  Terminal,  4010,  (2) 

Station  One  (Figure  No.  5-2) 

o  42"  X  60"  active  Area  Data  Automation  (X/Y  Digitizer  Table) 
o  Standard  Cursor  (five  push  buttons) 
o  Tektronix  4010  CRT 
o  16  Key  Keyboard 

Station  Two  (Figure  No.  5-3) 

o  42"  X  60"  active  area  Data  Automation  (X/Y  Digitizer  Table) 
o  Special  cursor  (LED  display/five  push  buttons) 
o  Tektronix  4010  CRT 

o  Threshold  Technology,  Inc.,  -  Model  500  voice  data  entty  temnn 
o  16  Key  Keyboard 


SERIAL  •  BIT  DATA  BSOO  BAUD 


Figure  No.  5-2 


5 .  3  BDRS  Software  Con  f  icjurat  ion 


5.3.1  Key  Areas 

The  key  areas  which  pertain  to  the  utilization  of  vendor  supplied 
system  software  (MRDOS/INFOS)  within  the  BDRS  are  as  follows: 

o  Multitasking  -  A  multiple  task  environment  is  one  in  which  logically 
distinct  tasks  compete  simultaneously  tor  the  use  of 
system  resources. 

o  Overlays  -  Overlaying  is  a  technique  used  for  loading  routines 

into  main  memory  from  some  type  of  mass  storage  during 
the  execution  of  a  program. 

o  User  Device  -  User  devices  are  any  device  that  is  not  part  of  Data 
General's  equipment  identification  during  a  system 
generation  (SYSGEN) .  These  devices  are  identified 
to  the  system  at  run  time  via  operating  system 
subprogram  calls.  (Reference  Data  General  Application 
Note:  "User  Device  Driver  Implementation  In  The  Real 
Time  Disk  Operating  System",  017000002-03). 

5. 3. 1.1  Multitasking 

The  BDRS  subsystem  software  utilizes  the  multitasking  capability  pro¬ 
vided  by  Data  General's  Mapped  Real  Time  Disc  Operating  System  (MRDOS/INFOS). 
The  digitization  software  consists  of  three  (3)  tasks  which  can  run  in  the 
foreground  of  the  C300  Data  General  processor.  The  first  task  functions 
merely  to  wake  up  the  two  remaining  tasks;  table  1  and  table  2.  Thus,  two 
digitizing  tables  can  be  operated  concurrently  in  one  ground  of  the  system. 

In  the  remaining  ground,  either  Data  Base  or  Batch  can  be  executed. 

5. 3. 1.2  Overlays 


The  table  1  and  table  2  software  modules  are  very  similar.  The  exe¬ 
cution  of  each  is  divided  into  overlays.  Each  task  has  one  overlay  segment 
assigned  to  it.  A  detailed  examination  of  the  actual  overlay  structure  by 


( 

overlay  iioinu  and  routine  name  is  depicted  in  Volume  I  Bathymetric  Data 
Reduction  Subsystem  Software  Documentation  BOC  Phase,  Section  3.2.2. 

Each  of  the  three  (3)  data  base  related  processing  modes  (On-Line,  Batch 
and  Master)  utilize  the  overlaying  technique  in  performing  tlieir  respective 
processes.  Consistent  througljout  the  above  mentioned  processes  is  the  fact 
that  each  processing  function  is  defined  as  an  overlay  and  is  loaded  into 
core  only  as  the  result  of  a  user  selecting  the  particular  function.  The 
implementation  of  such  a  philosophy  was  made  simple  by  the  fact  that  the 
tlirec  processing  modes  were  developed  in  a  modular  fasliion,  thus  allowing 
for  straight  forward  overlay  identification  and  structuring. 

5. 3. 1.3  User  Device  Implementation 

The  digitizing  tables,  voice  recognition  unit,  and  each  Telttronix 
4010  CRT  associated  with  a  table  are  introduced  to  the  system  software  in  ■ 

the  same  way.  The  method  of  introduction  is  described  in  detail  in  Data 

General's  User  Manual  entitled  "User  Device  Implementation".  It  suffices  i 

here  to  say  that  each  device  is  assigned  a  number  and  the  interrupt  service 
routine  address  for  each  device  is  stored  rn  a  system  vector  table. 

When  interrupt  occurs  in  a  multitasking  environment,  the  multi¬ 
tasking  activity  is  suspended.  When  the  interrupt  has  been  serviced,  it  is 
the  responsibility  of  the  interrupt  service  routine  to  reawaken  the  environ¬ 
ment.  This  is  accomplished  by  sending  the  appropriate  message  to  the  task 
interrupted  which  has  issued  a  receive  message  request  prior  to  the  interrupt. 

The  system  then  activates  this  task  and  the  multitasking  environment  is  , 

rescheduled. 

5.3.2  Conclusions  and  Recommendations 

Utilizing  the  INFOS/MRDOS  operating  system  (normally  called  INFOS)  both 
foreground  and  background  processing  may  be  operated  concurrently.  As 
stated  in  paragraph  5. 3. 1.1,  two  digitizing  tables  can  be  operated  in  one 
ground,  namely  the  foreground  and  either  the  Data  Base  or  Batch  functions 
can  be  operated  in  the  background.  The  operating  system  (INFOS),  digitizing 
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fulictioits  diiJ  tunL-t  ions,  an*  highly  ovorlayoii  piooossos,  that  is, 

they’  vt'ly  on  disk  I'osidont  ovorlays  for  loaviiinj  funotion.il  sul>pt osirams 
related  to  their  spooifio  fuiutions.  This,  ooupled  with  tl\o  faits  tluit 
the  ditjititin^  prooossos  is  a  time  oritioal  activity,  roiiuirinii  system 
resources  on  "demand"  (i.o.,  Ct'U)  while  data  collection  is  taking  place, 
and  the  BDRS  data  base  INfXtS  files  are  large  randomly  organised  disk  files 
which  need  be  accessed  an  indeterminate  number  of  times  depending  on  the 
complexity  of  the  user  requests  for  data}  supports  Synectic's  beliefs  that 
the  disk  is  over  burdened  due  to  the  frequency  and  location  of  disk  access 
which  caused  excess  head  movement  otherwise  known  as  thrashing. 

To  alleviate  the  above  mentioned  problems,  two  possible  solutioits 
are  recommet\dcd.  The  first  solution  would  be  to  dedicate  a  disk  (pre¬ 
ferably  a  quick  fixed  head  disk)  to  be  used  exclusively  for  storasje/ 
retrieval  of  the  following  data;  MRDOS/INPOS  overlays,  BDRS  subsystem's 
overlays,  intermediate  working  files  and  any  non-data  base  related  disk 
activities.  H\e  second  (Mssible  solution  would  be  to  dedicate  a  processor 
exclusively  for  the  data  base  fuitctions  sharing  a  master  disk  for  iitfre- 
quent  bulk  data  exchange  between  the  two  processors. 

It  should  also  be  pointed  out  that  in  the  production  envin.>iunent ,  the 
present  line  printer  output  cai'ability  may  be  inadequate.  TIte  Centronics 
Line  Printer  Model  101,  currently  part  of  the  BDRS  hardware  configuratioi^, 
outputs  165  characters  per  second  (cps) .  When  volume  hardcopy  reports  are 
necessary,  the  line  printer  and  disk  line  printer  spooler  are  tied  ui'  fox- 
long  periods  of  time.  When  this  occurs,  no  other  function  may  use  the 
printer  or  spooler  for  producing  any  hardcopy  output.  With  this  fact, 
serious  consideration  should  be  given  to  acquiring  a  faster  line  printer 
for  the  production  environment. 
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Rome  Air  Development  Center 

RAOC  pioM  and  executes  auiMch,  davctopment,  tzAt  and 
Aotectid  oc^tiuitLon.  pioffum  in  ^wppoAt  Command,  ContKoL 
CornmniAOtlonA  and  ItUelUgence  (C^j)  activUiu.  TcchnicaZ 
ana  enginteA^ng  AuppoKt  uiCtfUn  aA.ea6  o^  ttcJinicat  competence. 
4A  pfiovAded  to  ESO  Pnogaam  OUiceA  (PO4)  and  ESP 

tlenitn^.  The  pfUneipat  tecnnicat  mii6ion  oacm  a/te 
commnicationA,  etec^umagneXic  guidance  and  contAot,  40/1- 
vtUtance  of  gAound  and  aeAoopace  objeciA,  inteUigence  data 
cnltectioH  and  handtO.g,  infoAnation  igotem  teciinotogy, 
AonoAjpheAie.  pAopagation,  Aotid  itate  icienceA,  mccAouave 
pfufUcA  and  electAoncc  AeliablUtu,  maintainabilUu  and 
compatcbitUg. 


